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ABSTRACT 
Many  methods  have  been  proposed  and  tested  In  the 
past  to  obtain  an  optimal  response  for  step  Inputs  to 
automatic  feedback  control  systems.  Most  of  these 
methods  have  led  to  the  use  of  sophisticated  control 
devices  ranging  from  small  analog  to  large  digital 
computers.  Here  the  possibility  of  a  simplified 
switching  logic  combined  with  an  open-closed  loop 
servome onanism  is  investigated.   Deadbeat  response  to 
step  inputs  was  the  object  of  this  study  rather  than 
a  time  optimal  response.   Two  types  of  logic  were 
investigated.   A  time  invariant  controller  was  analyzed, 
built  and  tested.   The  system  works  on  the  principle 
of  constant  switching  times  with  the  output  being 
controlled  by  an  open  loop  driving  voltage  which  is 
proportional  to  the  input  step  size.  At  the  completion 
of  the  open  loop  mode  of  operations,  the  system  is 
returned  to  the  normal  closed  loop  mode. 

The  writers  wish  to  express  their  appreciation 
for  the  assistance  and  encouragement  given  them  by 
Dr.  John  Ward,  of  the  U.  S.  Naval  Postgraduate  Schools 
in  this  investigation. 
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1 .   Introduction. 

Numerous  schemes  have  been  applied  to  positioning 
feedback  control  mechanisms  to  obtain  optimum  response 
for  step  inputs.   Some  of  the  methods  used  have  been 
(1)  dual  mode  operation   using  a  relay  control  computer 

which  approximates  the  optimal  switching  llne9  (2)  use 

p 
of  compensation  networks  to  improve  relay  performance  , 

(3)  application  of  discontinuous  damping  to  a  relay 

servo^,  (4)  conditional  switching  techniques \    and 

others  5>6e   jn  general,  the  basic  idea  behind  these 

various  approaches  is  the  desire  to  obtain  deadbeat 

response  in  minimum  time  for  a  step  input. 

The  time  optimization  requirement  inherently 

requires  a  sophisticated  controller  that  can  compute 

the  proper  switching  times  accurately.  A  question 

#  *  *  *  -*  %  #  fl- 
IC.  C.  Matthews,  R.  C.  Boe,  The  application  of 
nonlinear  techniques  to  serv ©mechanisms,  National 
Electronics  Conference  Vol.  VIII,  pp.  10-21. 

2D.  McDonald,  Nonlinear  techniques  for  improving 
servo  performance,  National  Electronics  Conference, 
Vol.  VI,  pp.  400-421, 

■^Harris,  McDonald,  Thaler,  Quasi-optlmlzation  of 
relay  servos  by  use  of  discontinuous  damping.  Applica- 
tions and  Indus try 3   November,  1957. 

^S.  I.  Leberman,  A  bang-bang  attitude  control 
system  for  space  vehicles,  Aerospace  Engineering,  Oct.  1962 

-*T.  R.  Frederlckson,  A  time-optimal  positioning  servo, 
Control  Engineering,  February,  1963. 

G.  J.  Thaler,  M.  P.  Pastel,  Analysis  and  design  of 
Nonlinear  Feedback  Control  Systems.   Chapter  7. 


which  might  be  asked  is:   "If  the  requirements  for 

time  optimization  were  dropped,  retaining  the  provision 

for  deadbeat  response,  would  it  be  possible  to  design  a 

simple  controller  for  a  second  order  system  that  would 

be  of  practical  value?" 

In  this  paper  two  approaches  to  the  development  of 

a  simple  controller  for  a  second  order  system,  are 

investigated.   In  both  cases  the  system  operates  in 

two  modes,  linear  and  nonlinear,  which  correspond  to 

closed  and  open  loop  modes  respectively.   In  the  non~ 

linear  or  open  loop  mode  a  driving  voltage  Is  applied  to 

the  plant  for  a  prescribed  amount  of  time.   The  voltage 

is  then  reversed  until  the  system  reaches  an  output 

which  either  equals  or  approximates  the  input.  At  this 

point  the  controller  returns  the  system  to  the  linear, 

closed  loop,  operation. 

As  previously  stated,  in  both  approaches  invest!- 

gated,  relay  switching  is  caused  to  occur  on  a  time 

basis.   One  of  these  approaches  sets  the  switching 

times  proportional  to  the  command  signal  while  the  other 

holds  the  switching  times  constant  while  setting  the 

driving  voltage.  The  latter  scheme  is  suggested  by  the 

principle  of  superposition. 

The  system  of  Fig.  1  is  the  basic  block  diagram  for 

both  approaches.   The  step  input  commands  are  represented 

by  R,  the  system  output  position  by  C  and  the  output  of 

the  open  loop  controller  by  V.  The  plant  is  second  order 

„2- 


and  has  the  transfer  function: 


(1) 


K/f 


>(s+fr) 


The  system  Is  to  be  controlled  by  an  open  loop  controller 
which  functions  on  command  signals  only.  Any  load 
perturbations  or  other  similar  disturbances  within  the 
loop  will  not  affect  the  controller. 


OPE  A/    LOOP 
CONTROLLER 


Fig,  1  -  Basic  System 
The  assumptions  made  for  the  theoretical  investi- 
gations are  (1)  the  system  is  linear,  (2)  the  system  is 
second  order,  (3)  command  signals  are  to  be  steps  only 
and  (4)  new  commands  are  not  to  be  given  until  the 
previous  command,  has  been  completely  executed. 


2.   Time  proportional  to  command  signal. 

With  the  output  voltage  of  an  open  loop  controller 
a  constant  in  magnitude,  there  is  a  specific  relation- 
ship for  the  time  to  reverse  the  polarity  of  the  driving 
voltage  and  the  time  to  return  the  system  to  normal 
linear  operation  in  order  to  force  deadbeat  response. 
By  using  an  approximation  to  this  relationship,  a  theore- 
tical controller  was  derived.   The  system  Incorporating 
this  controller  was  then  investigated  on  the  CDC  digital 
computer. 

For  purposes  of  this  investigation,  the  proposed 
system  is  to  function  as  follows?   At  some  time,  t  ,  a. 
step  command  signal,  +  R,  is  applied  to  the  system. 
Application  of  the  command  signal  causes  the  loop  to 
be  opened  and  a  voltage,  +  V,  to  be  applied  to  the  plant 
for  a  period  of  time,  ts„  At  time  ts  the  polarity  of 
the  signal  is  reversed  for  a  period  of  time,  t  .  The 
total  time  of  nonlinear  operation,  tt,  is  defined  by; 

(2)     tt  =  tg  +  tp 
At  time  t,  the  controller  closes  the  loop  returning 
the  plant  to  standard  closed  loop  operation.   If  at 
time  t,  ,  the  plant  output  has  arrived  at  tbe  desired 
position  commanded  by  the  input,  then  true  deadbeat 
response  has  been  achieved. 

A  digital  program  was  written  to  compute  the  times 
ts  and  t^  for  various  values  of  inputs.   The  system 
parameters,  K  and  "T  ,  as  well  as  the  applied  voltage,  +V, 


were  held  constant.  For  the  system  under  evaluation 

the  following  values  were  arbitrarily  selected^ 

K/f    =  1,  1/V   =  5  and  V  =  100.   Each  combination 

of  tg  and  t^  that  results  in  deadbeat  response  specifies 

a  command  signal  R„  Thus^  from  the  program  data9  a 

plot  of  t  and  t,  versus  R  was  made  for  Q^R^TT        radians 

as  shown  in  Fig.  2,   From  these  data9  t  and  t  can  be 

s      t 

approximated  by  a  straight  line  as  shown  in  the  figure. 
These  straight  line  approximations  of  t  and  t.  can  be 
evaluated  in  point  slope  form  ass 

(3)  i>  =  .W|*|  +  0.05 

(4)  tt  '-  .092|>?l  +0J0S 

Since  ts  and  t^  are  independent  of  the  sign  of  the 
input,  the  magnitude  of  R  is  specified  in  the  above 
equations. 

It  is  immediately  obvious  from  Fig0  2  that  the 
straight  line  approximation  is  extremely  poor  for  small 
inputs,  but  this  problem  can  be  eliminated  if  the 
system  is  held  in  closed  loop  for  small  inputs.  Thus, 
in  this  system,  the  minimum  step  size  for  open  loop 
operation  would  be  about  ,25  radians 0 

Using  the  switching  times  determined  from  equations 
(3)  and  (4),  a  program  was  written  and  appended  to  "Program 
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controller.   Input  equals  1.5  radians. 
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Analog1   to  study  the  performance  of  a  system  utilizing 
this  type  controller*   Pour  runs  were  made  with  inputs 
of  «5»  1»  1.5*  and  2  radians.   The  resulting  phase 
planes  and  time  response  characteristics  as  plotted  by 
the  computer  are  shown  in  Figures  3  to  10,   By  Inspection 
of  Pig.  2  it  is  possible  to  predetermine  whether 
switching  will  occur  early  or  late*   For  example 9    if 
the  input  is  .5  radians ,  Pig.  2  predicts  that  both 
switching  times^  ts  and  tr9   will  be  late  and  the  system 
will  overshoot*  This  is  verified  in  Figures  3  and  40 
If  the  input  is  1.5  radians,  it  would  be  predicted 
that  switching  would  be  early.  This  is  verified  in 
Figures  7  and  8. 

This  approach  does  offer  some  Improvement  in 
time-to-steady-state  over  normal  linear  operation  but 
it  is  limited  to  a  comparatively  small  range  of  Inputs. 
A  controller  that  could  vary  the  times  of  switching  as 
specified  by  equations  (3)  and  (4)  would  not  be  simple 
to  constructs  which  is  contrary  to  one  of  the  main 
design  objectives  as  set  forth  in  section  1. 

'"Program  Analog"  was  written  by  Dr.  J»  Re  Ward 
of  the  U.  S.  Naval  Postgraduate  School.   The  program  as 
used  is  shown  in  Appendix  I. 


3.   Time  invariant  controller. 

Assume  that  a  given  linear  second  order  system  is 
controlled  in  such  a  manner  as  to  provide  a.  deadbeat 
response  for  a  step  input.   If  the  output  of  the  linear 
system  is  proportional  to  a  step  inputs  then  by  the 

o 

principle  of  superposition  9    the  output  is  doubled 
when  the  input  is  doubled.  Furthermore 9    if  the 
acceleration  time,  ts,  and  the  deceleration  time  t  , 
are  known  for  a  specified  step  commands  then  it 
will  be  possible  to  obtain  deadbeat  response  for  all 
step  commands  if  the  switching  times  are  held  constant 
while  the  driving  voltage  is  varied  in  proportion  to 
the  step  command. 

To  determine  the  appropriate  switching  times 0   it 
is  first  noted  that  if  the  controller  voltage  is  V9 
then  the  system  output  willg  if  the  initial  conditions 
are  zero,  be^ 

(5)    Cm-     KV/« 


S2($-*y<r) 


The  step  voltage  +V  is  to  be  applied  for  a  finite  time* 
tggl  and  the  step  voltage  -V  is  applied  for  the  finite 
time,  tr.   If  at  tg  +  t  ,  the  output  of  the  system  is 
identically  equal  to  the  input s    and  the  output  rate  is 
zero,  then  system  response  is  indeed  deadbeat.   It  is 

&  #  «  #  #  ■»{■■&  * 

o 

M.  E.  Van  Valkenburgg  Network  Analysis 9  Prentice- 
Hall,  Inc.,  pp.  79=80. 
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now  necessary  to  find  the  required  conditions  for  such 
deadbeat  responseo 

The  inverse  transform  of  equation  (5)  is; 

(6)    a**   =  M^t       +  U   -  <) 


To  obtain  the  system  rate  at  time  t  it  is  only  necessary 


to  take  the  derivative  of  equation  (6).   Thuss 


(7)      CM    *   KV(l  -€m¥r) 


The  reversed  voltage  Is  now  applied  for  a  period  tr9 
and  considering  the  Initial  conditions  that  exist  at 
time  t  9    it  follows  that  the  equations  for  system 
position  and  rate  for  ts  £  t£  t^  ares 

(8)   G>)-       vkA    +  _£^L.  +     cm  <W< 

SVS  +  '/f}       (S  +  *A)         S(S  +  l/f)         S(S+ ft) 


and; 


(9)    Cfs\z         W/f      ,   ,._gfrJL_ 
S(5  +  W       (3  +  Vt) 


ors 


»)    s   MCI*,  -*r  +<[/  -^        f2^        'J/ 


(10)   Cfc 


and; 

do  <r^)  s  wl-J +  £      (2 -*     U 

For  deadbeat  response  the  system  rate?  C  , 
must  equal  zero  at  time  t+„  Therefore?,  by  setting 
equation  (11)  to  zero,  the  exact  relationship  between  the 


acceleration  time,  t  ,  and  deceleration  tlrae»  t  9    is 


obtained j  namely? 

9  -  £ 


(12)    6     =2 


To  obtain  the  deceleration  time  required  for  any 
acceleration  time?  It  is  only  necessary  to  take  the 
natural  logarithm  of  both  sides  of  equation  (12); 

(12a)    ir    ~     fJ*(Z~C         > 

To  obtain  the  system  output  for  deadbeat  response 
in  terms  of  tq  and  tr,  equation  (112)  is  substituted 
into  equation  (10)  which  yields? 

(13)    CttA     =  VKftt-tr) 

Finally 9  the  system  output  in  terms  of  the  acceleration 
time  only  can  be  found  by  substituting  equation  (12a) 
into  (13)  to  give? 


(14)   CM    s 


„• 


is  -  fA 


Thus  by  using  a  controller  voltage  proportional 
to  R,  and  by  choosing  ts  and  tr  as  calculated  from 
equations  (12a)  and  (14)9  deadbeat  response  can  be 
achieved  for  any  step  input  provided  that  the  plant 

is  adequately  described  by  a  (linear)  second  order 
transfer  function. 
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4.   A  design  example  of  the  time  Invariant  controller* 

A  logical  sequence  of  manipulations  to  determine 
the  proper  switching  times  for  a  given  linear  second 
order  system  can  be  easily  obtained  from  equations 
(12a)  and  (14).   For  purposes  of  illustration 9    assume 
that  the  system  has  the  following  characteristics? 

Ks  10     Maximum  output  desired  =  10 

<f  ~  1»0    Saturation  will  occur  If  K?  is  greater 

than  100* 

Using  equation  (14)  with  the  constants  given  above 

and  selecting  a  maximum  applied  controller  voltage 
of  10  (to  avoid  the  saturation  non-linearity)  a  trans~ 
cendental  equation  for  the  acceleration  time^  tgs  is 
obtained? 


(14)   iO   »  /Oo[t*    -/•.(2-e*"*)J 


A  value  of  tg  equal  to  9368  seconds  will  satisfy  the 
equation*   It  is  then  necessary  to  find  the  deceleration 
time  tr.  To  do  this9  substltue  time  tg  into  equation  (12a) 

MB' 


(12a)  tr   ="  fn(  Z  -  £~     '  )  s  >26S    see, 

The  total  time  of  open  loop  operation  is  then  e636  seconds, 

Using  the  above  constants 9   the  system  equations  were 
solved  on  the  digital  computer   Only  one  input  was 
given  to  the  system  as  other  inputs  would  provide  the 
same  type  of  output  changed  by  a  constant  only8  Fig,  11 
shows  the  computed  phase  plane  response  of  the  system* 
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5.,   Test  of  the  time  invariant  controller  applied  to 
a  d-c  motor  plant s 

The  controller  analyzed  in  sections  3  and  4  above, 
is  here  Investigated,  experimentally  in  relation  to  a 
DC  motor  plant  which  was  predominantly  second  order 
and  linear  up  to  the  saturation  limit  of  the  driving 
amplifier*  Runs  were  made  with  various  step  inputs  in 
both  the  tinie= invariant  and  standard  closed  loop  modes. 

The  following  operations  by  the  open  loop  controller 
are  required  !   4  st  the  theory: 

1 )  At  time  tQ  a  step  input  is  eommandede  The 
loop  is  opened  and  an  open  loop  driving  voltage  is 
applied  to  the  plant , 

2)  At  time  tg  the  driving  voltage  is  reversed* 

3)  At  time  t^  the  driving  voltage  is  removed 
from  the  plant 9    the  loop  is  closed s   and  the  system 
returned  to  its  standard  closed  loop  mode  of  operation. 
The  operations  required  were  carried  out  by  a  timing 
device  consisting  of  a  variable  speed  D©  Co  motor 
geared  to  a  shaft  containing  three  8'make  and  break" 
switches.  The  time  of  contact  of  these  switches  rela<= 
tlve  to  each  other  was  variable^,  that  is9  their  relative 
anglular  positions  on  the  shaft  could  be  varied*  The 
switches  were  connected  to  relay  circuits  which  accomplished 
the  desired  switchings,  Thus9  relay  switching  times 

could  be  varied  hy   adjusting  the  speed  of  the  Do  G. 
motor  or  by  adjusting  the  relative  angular  position  of 
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the  "make  and  break00  contacts  • 

Fig*  12  shows  a  block  diagram  of  the  experimental 
D,  C8  servo  used  for  testing*   When  a  step  command  is 
received ,  relay  one  switches  to  its  normally  open 
position.   This  is  time  t0.  Relay  two  switches  to  its 
normally  open  position  at  time  tg  to  reverse  the  open 
loop  driving  voltage •  A  third  relay g   not  shown  in  the 
figure j,  is  used  to  return  relays  one  and  two  to  their 
normally  closed  positions  at  time  t+. 
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Fig.  12  D»  Go  Servo  with  a  time  Invariant  controller 
The  transfer  function  of  the  open  loop  system  was 
evaluated  ass 


(17) 


Z   Mr/ 1.8 


SfS  *  i/i.») 

where  the  gain  is  2«94  and  the  system  time  constant  is  1.8 

seconds  *  The  system  output  position^  C9   was  obtained  by 

integrating  the  output  of  a  tachometer  geared  to  the 
driving  motor.  As  a  result 9    the  system  output  was 
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evaluated  In  volts 9  and  consequently 9  the  command  and 

output  signals  could  be  compared  In  volts  rather  than 
radians*  Further  Investigation  of  the  open  loop 

system  revealed  that  an  open  loop  command  signal  of 
about  eight  volts  would  cause  saturation  In  the  transistor 
power  amplifiere  A  schematic  of  this  amplifier  is  shown 
In  Fig,  13. 

For  the  purpose  of  steady  state  accuracy 9  of 
primary  importance  in  positioning  systems^  the  closed 
loop  gain  was  set  at  29 «4  or  ten  times  greater  than  the 
open  loop  gain* 

To  simplify  the  construction  of  the  controller 9  it 
was  decided  to  obtain  switching  times  such  that  a  one 
volt  Input 9  R9   would  cause  a  one  volt  driving  voltage 9   V. 
In  other  words 9   the  constant  of  proportionality  derived 
in  equation  (14)  was  made  equal  to  one8  This  means  that 
in  equation  (14)9  V  =  C9  and  tg  can  be  evaluated 
directly  from  that  equation  to  yield  the  desired  value* 

(14)    Cut)    -   //cfc-//.fW"**)] 

J_  -  tt   -  /.8A(2-£    ) 

6  S  M  sec. 

By  using  the  value  of  tg  equal  to  «98  sec*,  In  equation 
(12a) 9  t  is  found  to  be  064  sec*  The  relay  switching 
times  were  then  set  accordingly* 

An  input  of  three  volts  was  used  to  check  out  the 
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switching  time  computed  aboveo  The  response  was  not 
deadbeato  The  switching  times  were  then  adjusted  to 
obtain  deadbeat  response  t©  the  three  volt  input©  The 
actual  time  of  voltage  reversal  was  11  a  11 5  see*  compared 
to  the  computed  e98  sec0  The  total  time  of  operation 
was  1W65  sec*  compared  to  the  computed  time  of  11 0  52  see* 
These  differences  were  probably  due  primarily  to  the  fact 
that  the  system  was  ©Ely  quasi  linear  and  second  order 9 
and  the  sj  .  m  int  were  determined 

assuming  a  linear  second  order  systems 

A  six  channel  Brush  recorder  was  used  to  record 
the  closed  loop  command  signal 9  the  ©pen  loop  command 
signal 9   the  system  output  and  the  system  rate0  Runs 
were  made  with  the  inputs  varying  from  1,4  volts  to  10 
volts*  A  standard  closed  loop  response  to  the  same 
step  commands  was  also  made  for  comparison  purposes • 
Figures  14  through  21  show  the  responses  obtained©  Each 
page  contains  a  tlme=  invariant  system  response  followed 
immediately  by  a  standard  closed  loop  response  to  the 
same  signal©  Table  1  is  a  comparison  of  the  different 
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system  with  a   1.94  volt  input. 
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"Fig.  16.  Time  invariant  and  closed  loop  response  of  DC  servo 
system  with  a  4.0  volt  input* 
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The  controller  reversed  the  input  to  the  motor  at  the 
proper  time  but  the  reversal  is  not  indicated  on  this 
trace.  It  was  not  possible  to  record  the  input  to 
the  motor  (amplifier  output)  showing  the  reversal  of 
driving  voltage*  because  of  grounding  problems0   A 
ground  at  the  motor  input  would  cause  one  or  both  of 
the  power  transistors  to  burn  out  in  the  power  ampli- 
fier. The  third  trace  is  the  system  output  and  the 
fourth  trace s  the  system  rate.  The  time  ts,  the  time 
at  which  the  polarity  or  sign  of  V  is  reversed^  is 
readily  observed  on  the  last  trace,  system  rate9  as 
there  is  a  definite  change  in  slope  at  that  time* 

Time  on  all  the  traces  reads  from  left  to  right  with  t 

°        o 

signified  by  the  coanencement  of  the  open  loop  command 
signal  and  t.  by  the  end  of  this  command  and  the  start 
of  the  closed  loop  command  signal. 

It  was  mentioned  earlier  that  an  output  signal  of 
eight  volts  would  saturate  the  power  amplifier*.   It 
follows,  then,  that  if   C>^   ,  the  amplifier  will 
saturate  when  the  driving  voltage  is  reversed „  The 
reason  for  this  is  that  at  the  time  of  voltage  reversal 
a  back  electromotive  force  has  built  up  in  the  motor 
which  approaches  the  driving  voltage.  This  back  emf 
has  the  effect  of  increasing  current  requirement  from 
the  amplifier  when  the  driving  voltage  is  reversed. 
Therefore,  at  the  time  of  switching,  the  system  should 
tend  toward  saturation  if  the  signal  voltage  is  greater 
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than  four  volts  *  This  saturation  caused  a  reduction 
in  the  braking  or  deceleration  power  available 0  From 
this  it  follows  that  the  system  output  would  tend  to 
overshoot.  These  effects  may  be  seen  in  Figures  116, 
17,  and  18,  By   analyzing  the  system  rates  at  time  tg 
and  t+  for  these  figures  it  can  be  demonstrated  that 
saturation  does  occur  in  the  time  increment  ts  S  t  is,   t^. 
At  time  ts,  for  the  three  figures  mentioned  above s    the 
system  rate  increased  proportionately  with  the  command 
signal  as  called  for  by  the  theory.  Yet  at  time  t^9 
the  system  rate  was  larger  than  predicted  meaning  that 
the  system  did  not  decelerate  sufficiently  which 
resulted  in  overshoot  at  time  t. .  For  example 9    the 
system  rates  at  time  t  and  t.  for  the  five  volt  Input 
(Fig,  17)  are  6,1  and  ,95  volts  per  second  respectively , 
For  an  input  of  6,2  volts  (Fig,  18)  one  would  predict 
by  the  principle  of  superposition  that  the  rates  would 
be  increased  by  the  factor  6,2/5.0  or  would  be  7o5 
and  1.1  volts  per  second.  Actually »  Fig,  18  shows 
the  rates  to  be  7»5  and  1,5  volts  per  second.  The 
rate  at  time  tg  agrees  with  the  theory  but  at  time  t^ 
the  rate  differs  by  *4  volts  per  second.  This  difference 
can  be  accounted  for  by  the  saturation  that  occurs 
when  the  voltage  Is  reversed. 

When  the  system  input  exceeds  eight  volts 9    saturation 
is  predicted  at  t0  and  at  ts  ,   One  would  expect 9  there° 
fore,  that  the  system  output  would  never  reach  the 
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desired  point  at  time  t^  or  in  other  words ,  the  system 
will  have  undershoot  at  time  tt.  The  reason  for  this 
phenomenon  is  that  the  driving  voltage  is  never  large 
enough  to  provide  the  response  called  for  by  the 
theoretical  equations.  The  results  of  this  type  of 
saturation  are  shown  in  Figures  19  and  20.  Fig.  20 
clearly  demonstrates  the  effect  on  system  output. 

Pig.  21  shows  the  response  of  the  system  with  a 
small  input,  1.4  volts.  Here  the  system  output^  by 
time  t^,  does  not  reach  the  desired  value0  The  primary 
cause  for  this  is  that  the  small  driving  voltage  is  of 
insufficient  magnitude  to  overcome  the  stiction  and 
friction  in  the  system* 

Prom  these  experimental  tests  it  can  be  seen  that 
the  time  invariant  switching  scheme  does  provide 
deadbeat  response  over  a  range  of  inputs.  If  deadbeat 
response  is  desired  in  a  positioning  feedback  control 
system,  the  method  of  time  invariant  switching  can 
provide  this  in  a  given  time  increment  over  a  designed 
range  of  inputs. 
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6.   Test  of  the  time  Invariant  controller  applied  to 
an  amplldyne  driven  motor. 

The  time  invariant  controller  has  been  shown  to 
work  over  a  range  of  Inputs  with  a  quasi  linear 9 
second  order  system.  One  might  now  inquire  as  to  its 
usefulness  with  a  system  that  Is  quite  nonlinear  in  that 
it  exhibits  hysteresis^  excessive  stiction9  and  satura- 
tion.  To  investigate  this,  a  system  was  constructed 
that  incorporated  an  amplldyne,  a  £  HP  d-c  shunt  wound 
motor,  and  a  permanent  magnet.  Do  0,  generator  The 
output  of  the  generators  or  tachometer^  was  integrated 
to  simulate  the  system  output.  As  in  the  previous 
experiment  the  system  output  and  input  was  in  volts „ 
Pig.  22  depicts  the  system  as  set  up  in  the  laboratory 
with  the  relays  performing  the  same  functions  as  out- 
lined in  section  5,  Fig,  23  is  a  picture  of  the  system 
as  it  appeared  in  the  laboratory. 


open  loop 

CONTROLLER 


Pig,  22  Amplldyne  driven  motor  servo  system  with 
time  invariant  controller. 
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A  plot  of  the  voltage  input  versus  the  tachometer 
voltage  output  for  the  open  loop  system  appears  on 
Fig.  24.   Notice  that  at  least  two  volts  are  required 
to  start  the  system  in  motion  and  that  at  about  six 
volts  the  system  starts  to  go  into  saturation,.  The 
hysteresis  loop  is  rather  broad  which  could  provide 
difficulty  in  system  performance.  Prom  open  and  closed 
loop  tests  of  the  system,  without  the  time  invariant 
switching,  the  open  loop  transfer  function  was  determined 
to  be  approximately s 

(17)  -iQ 

S(S  +  7) 

Depending  on  the  test  being  conducted*  the  system  gain 

varied  from  100  to  115  and  the  reciprocal  time  constant 

from  6.9  to  7.5  inverse  seconds. 

Using  the  values  of  gain  and  the  time  constant 

from  equation  07)  and  having  a  unity  relationship  between 

the  open  loop  output  and  open  loop  driving  vcltage?  t_ 

was  determined  from  equation  (14)  to  be  about  033 

seconds.  Then  t  9  from  equation  (12a)9  was  evaluated 

as  .27  seconds.  The  system  was  run  with  an  input  step 

of  four  volts.  The  response  was  not  deadbeat  so  the 

times  of  switching  were  adjusted  to  yield  a  near  dead= 

beat  response.  The  timing  circuit  controls  were  not 

fine  enough  to  permit  precise  adjustments  of  the  switching 

times.  The  adjusted  times  were  then  measured  and  found 

to  be  .38  second  for  time  tg  and  .09  second  for  tr. 
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Fig.    24.        Hysteresis    loop   of  amplidyne   driven  servo   system. 
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The  acceleration  time  was  ,05  second  longer  than  the 
computed  value  and  the  deceleration  time  was  .18  second 
shorter  than  the  computed  value.  This  shorter  deeelera- 
tion  time  is  probably  due  to  the  large  friction  that 
was  present  in  the  motor* 

A   Brush  recorder  was  used  to  record  the  values  of 
the  closed  loop  command  signal,  the  open  loop  command 
signal,  the  system  output,  the  system  rate^,  and  the 
voltage  applied  to  the  motor.  These  recordings  are 
shown  in  Figures  25  through  38«  Following  each  time=> 
invariant  run  to  a  given  input,  the  system  was  run 
in  the  normal  closed  loop  mode  with  the  same  quantities 
being  recorded.  A  figure  of  the  closed  loop  recording 
follows  each  time  invariant  figure.  The  maximum  over= 
shoot  and  time~to=>steady*=state  was  evaluated  for  all 
runs  and  a  comparison  of  the  times  to  steady  state  was 
made.   The  results  of  this  evaluation  are  tabulated  in 
Table  II. 
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TABLE  II 
.£***************************.***   # 

Amplidyne        Open-closed  loop  *       Closed  loop 

*  Servo  *  * 
****************************** 

*  Input  Mpt'      t  *   *  Mpt     tgs   At*      * 

********************&*****  *°*  *  * 
*  * 

3.0  volts        1o0     1.18  sec.   1.0   059  sec.  <=.59  sec. 

*  * 

3.5  KO     1,00       1eO   .50      °.50 

*  * 

4.0  KO      .50        1.26  1.00       .50 

4.5  1.0      .48       1.31  1.05       .57 

*  * 

5.0  1.0      .48       1.58  1.00       .52 

5.5  U25    1.03        2.02   e93      -.10 

*  * 

6.0  1.88    1.00        2o12  1.30       .3 

********************************* 


Note; 

1.  M^s   Ratio  of  maximum  output  to  desired  output, 

2.  Time  to  steady  state  operation^,  taken  to  be 
when  system  is  within  ,1  volts  of  command 
signal. 

3.  At-ss  is  defined  as  the  difference  between 
the  times  to  steady  state  for  the  open-elosed 
loop  and  closed  loop  operation. 

Prom  the  recordings  it  can  be  seen  that  the  range 

from  four  to  five  volts  the  response  was  almost  dead~ 

beat.  (Figures  25  through  30.)   When  the  input  voltage 

was  in  excess  of  five  volts  the  system  output  started 

to  overshoot  the  desired  value  at  time  t^.  As  in 
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Fig.  25.   Time 'invariant  response  of  amplidyne  driven  motor  servo 
system  with  a  four  volt  input* 
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Fig*  26«   Closed  loop  response  of.  amplidyne  driven  motor  servo 
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Fig.  27.   Time  invariant  response  of  amplidyne  driven  motor  servo 
system  with  a  4|-  volt  input. 
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Fig.  28.   Closed  loop  response  of  amplidyne  driven  motor  servo 
system  with  a  -4^-.  volt- input. 
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Fig.  29«   Time  invariant  response  of  amplidyne  driven  motor  servo 
-system  with- a-  five -volt- -input  •-■ 
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section  5  with  the  Do  Co  servo 9    this  is  due  to  satura- 
tion of  the  amplidyne  when  the  driving  voltage  is 
reversed.  This  type  of  response  is  shown  in  Figures  31 
through  34. 

For  step  inputs  less  than  four  volts  system 
performance  deteriorated  rapidly.  This  was  due  almost 
entirely  to  the  excessive  stiction  and  friction  in  the 
system.   It  was  found  that  the  system  did  not  respond  to 
a  command  less  than  two  volts  or  about  30%   of  the 
operating  range.  The  system  response  to  inputs  of  3 
and  3»5  volts  are  shown  in  Figures  35  through  38* 

As  might  be  expected ,  the  use  of  an  open-loop 
controller  with  a  far»from° ideal  plant  does  not  lead 
to  acceptable  performance. 


5q 


~Time~« 


71        72~T5       74     ~  5: 


1.0      ceconds 


VOLTS 


6 


2 


VOLTS 


Closed  loop  command  signal 


6 
4 
2 

-0- 


CO 


Open  loop  command  signal 


Input  to  d-c  motor 


Fig.  >1.   Time  invariant  response  of  amplidyne  driven  motor  servo 
system  with  a  5f  volt  input.-- 
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Fig.  52.   Closed  loop  response  of  amplidyne  driven  motor  servo 
system  with -a  5^-  volt-  input. 
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Fig*    3>      Time   invariant  response   of  amplidyne   driven  servo 
system  with -a -six  volt    input* 
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Fig.  54.   Closed  loop  response  of  amplidyne  driven  servo 
..system -with-  a  six  volt  input.  - 
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Fig*   55*      Time   invariant  response   of  amplidyne   driven  motor   servo 
-system-with  three  volts-  input. 
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Fig.    56.      Closed   loop  response    of   ampiidyne   driven  motor    servo 
systeauwith-a-three  volt   input.- 
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Fig*  57  •   Time  invariant  response  to  amplidyne  driven  motor  servo 
system  with  a  5s  volt  input* 
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7.   Conclusions. 

In  the  search  for  a  simple  scheme  to  provide 
deadbeat  response  for  step  Inputs  to  a  second  order 
system,  attention  was  concentrated  on  an  open-closed 
loop  controller.   For  the  system  investigated*  switching 
time  in  the  open  loop  mode  were  fixed  and  the  output 
from  the  controller  to  the  plant  was  proportional  to 
the  magnitude  of  the  step  input. 

Using  this  scheme,  any  linear  system  can  be 
controlled  with  deadbeat  response  (Principle  of  super- 
position).  To  control  the  second  order  plant  considered 
in  this  study  only  one  switching  point  is  required. 
The  experimental  program  was  intended  to  suggest  the 
shortcomings  of  this  approach  when  applied  to  a  relay 
system  with  its  non-linearities  and  higher  order  dynamics. 

The  first  real  system  studied  was  relatively  free 

from  any  non-linearities  and  provided  near  deadbeat 

response  over  approximately  50$  of  the  designed  operating 

range.  A  system  with  rather  pronounced  non-linearities 

was  also  chosen  for  study  and  provided  near  deadbeat 

response  over  less  than  20$  of  the  desired  operating 

range.   It  is  concluded ,  then,  that  if  deadbeat  or  near 

deadbeat  response  is  desired,  the  time  invariant  controller 

investigated  in  this  thesis  can  be  used  to  force  this 

response  in  a  linear  system  •   With  a  non-linear  system, 

the  controller  is  of  little  use.   In  any  event ,  the  time 

of  response  to  any  size  step  requires  the  same  Interval 

of  time  for  open  loop  operation. 
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APPENDIX  I 
PROGRAM  ANALOG 
A  Control  Data  Corporation  CDC  11604  digital 
computer  is  available  at  the  Ue  S.  Naval  Postgraduate 
School.   Dr.  J.  R.  Ward  of  the  Postgraduate  School^ 
Department  of  Electrical  Engineering 9  has  written  a 

Q 

FORTRAN  language  program  for  this  computer  which 
simulates  an  analog  computer.  The  program  is  referred 
to  as  Program  Analog.   It  retains  all  the  versatility 
of  the  analog  computer  and  provides  both  a  graphical 
and  numerical  output.  With  this  program  many  of  the 
problems  of  magnitude  and  time  sealing  as  well  as  the 
difficulties  In  simulating  nonlinearitles  associated 
with  an  analog  computer  are  avoided.  The  program j,  which 
is  reproduced  in  Figures  39  through  41 g   uses  a  modified 
Runge-Kutta  integration  scheme* 

This  program  was  used  for  the  investigations 
carried  out  in  section  2e   Fig.  39  Is  the  main  part  of 
the  program.   It  consists  of  the  necessary  instructions 
for  the  analog  simulation  of  a  given  problem  and  does 
not  require  any  changes  as  the  problem  is  varied .   This 
part  of  the  program  may  be  thought  of  as  the  analog 
computer.   Fig.  40  Is  the  FORTRAN  simulation  of  the 

Ik        5v   tx        ff^r   *t   r&        &v        o& 

°D.  D.  McCracken9  A  guide  In  FORTRAN  programming s 
John  Wiley  &  Sons9  Inc.,  New  York9  London9  1961, 
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problem  board.  This  part  of  the  program  must  be 
changed  in  the  same  manner  as  one  would  change  the 
wiring  on  a  problem  board*  Fig.,  411 9   the  data  Input9 
sets  the  initial  conditions 9   potentiometer  settings,, 
problem  time  and  provided  means  by  which  various  output 
data  can  be  obtained,  A  typical  numerical  output 
for  a  Program  Analog  problem  appears  in  Fig*  42 « 
Figures  3  through  10  in  the  body  of  the  th  -  are 
typical  outputs  for  the  program* 
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Fig*  59  •      Program  Analog. 


.JOB    MALLEY     (     10    NINITES    SHOULD    CO) 
PROGRAM    ANALCG 
CALL    EXEC 


END 

SUBROUTINE  EXEC 

DIMENSION  X(15C),XCOT(10  0),C(15C),ITITLE(10),JTITLE(10),KTITLE(4G) 

1  ,IT(5),IPn0),IG(10),CT(5)tXO(100),PR(10),CR(lC),Xl  (900)  , 

2  Yl(90C),X2(9  00),Y2(90  0),X3(9C0),Y3(900),X4(9CC) , Y4(900), 

3  X5(90C) ,Y5(900 ) ,A( 100) 
DO  1  J=  1  , 1 50 
X{ J)  =  C. 


1  C(J)  =  0. 
DO  1  CO  1  J  =  1,10C 

iooi  xo( j  )  =  o. 

DO  1002  J=l,  1C 
ITITLE! J)  =  8H 

1002  JTITLE(J)  =  EH 
DO  1003  J=l  ,  40 

1003  KTITLE(J)  =  8H 
NRC  =  0 
PRINT  200 

200  FORMAT! 1H1 ) 

201  FORMAT!/) 

202  FORMAT!//) 

203  FORMAT! ///) 
REAC  102,  !  ITITLE!  J),J  =  1  ,5) 

102  FORMAT! 10A8) 

READ  101, N, ITEST 
101  FORMAT!  12, A8) 

ICHECK  =  8H  EQUATIC 

IF! ICHECK  -  ITEST)  2,3,2 

2  PRINT  204 

204  FORMAT! 18H  DATA  FORMAT  ERROR) 
STOP  1 

3  READ  101 ,NR, ITEST 
ICHECK  =  8H  RUNS 
IF! ICHECK  -  ITEST)4,6,4 

4  ICHECK  =  8H  RUN 
IF! ICHECK  -  ITEST)  7,6,7 

6  IF(NR  -  9)  3,8,7 

7  PRINT  2C4 
STOP  2 

8  READ  101, NT,  ITEST 
ICHECK  =  8H  TITLE  C 
IF! ICHECK  -  ITEST)  10,9, 10 

9  IF(NT  -  9)11,11,10 

10  PRINT  204 
STOP  3 

11  DO  12  1=1, NT 
READ  ' 

12  PRINT 
IF!NR 

13  PRINT  206 
GO  TC  15 

14  PRINT  207, NR 

205  F0RMAT(5X,  1CA3) 

206  FORMAT! 5X,20hl  RUN  IS  CALLED  FOR.) 

207  FORMAT! 5X, II  ,2  IK  RUNS  ARE  CALLED  FOR, 

15  PRINT  203 
300  NRC  =  NRC  +  1 

READ  102, ITEST 

ICHECK1  =  8HZFRC  CCE 

ICHECK2  =  8HFCLC  CCE 

IF! ICHECK1  -  ITEST)  16, 18 , 16 

16  IFUCHECK2  -  ITEST)  17,20  ,17 

17  PRINT  204 
STOP  4 

18  DO  19  J=l , IOC 

19  C(J)  =  0. 

20  REAC  101  ,NC,  ITEST 

=  63- 


T  _  1    MX 

102, (JTITLE! J), J=l, 10) 

205, (JTITLE! J) ,J  =  1 , 1C) 
.  -  1)10,13,14 

206 
!  15 
PRINT  207, NR 
F0RMAT(5X,  1CA3) 
FORMAT! 5X,20F1  RUN  I! 
FORMAT( 5X, II ,21K  RUN: 
PRINT  203 


Fig.  59  .  Cont. 
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22 
23 

103 


ICHECK  =  8H 

IF(  ICHECK  - 
PRIM  2CL 
STOP 


CfEFF. 
ITEST) 


2  1,22,21 


STOP  b 

IF(NC)2  3,2  5,2  3 

UO  2L  J=1,NC 

REAC  103,  (  IT  (K),CK 

FORMAT (  5(IU,F12.U)) 

DO  2U  K"l ,5 

ITK  =  IKK) 


K)  ,K=  1  ,5) 


,  (K) 

"  FCrKTAT(12HNRLN  NUMBER  ,I1,///,U6H  THE  NON-ZERO  DATA  COEFFICIENTS  A 
IRE  AS  FOLLOWS,/) 

=     0' 


25 
208 


26 
209 


K    =    0' 
DC    27    J=l ,99 
IFIC  (J)  )26, 27,26 
PRINT    209, J,C( J) 
FCRVAK  1CX2H.C(  ,  12, 
K    =    K    +     1 

^AkiT  TklllT 


UH)  =  ,E12.5) 


27  CONTINUE         ,^rt_ 
IF(K)  1027,  1026, 1027 

1026  PRINT  240 

2U0    FORMAK  5XiiH\CNF^ 

1027  PRINT  202 

r,  r-    A  r^  1  O  ")    TTCCT 


READ  102, ITEST 
ICHECK1  =  8HZERC 
ICHECK2  =  8HFCLC 
ICFT 


ICS 
ICS 


ICHECK^  =  «Hhl_LL   1L5 

ICHECK3  =  8HFCLC  FCS 

IF( ICHECK1  -  ITEST)iC,28,3C 

DO  29  J=1,N 

XC( J )  =  0. 

IF{ICFlicK2  -  ITEST)31,35,31 

IFUCHECK3  -  ITlST)  i2,3i,3^ 

PRINT  2CU 

STOP  6 

DO  3U  J  =  1,N 

XC( J)  =  X( J)   _ 

READ  101, NI, ITEST 

ICHECK  =  3H  IC  CARD 

IFUCHECK  -  ITEST)36,37,  36 

PRINT  2CU 

STOP  7 

IF(NI)38,UC,38 

R°AD9103!uT(K),CT(K),K=l,5) 

DC  39  K=l ,5 
ITK  =  IKK) 
X0( ITK)  =  CT(K) 

F0RMAT<U8H  ThE  NCN-ZERO  INITITAL  CONDITIONS  ARE  AS  FOLLOWS,/) 

K  =  0 

DO  U2  J=l ,N 

IF(XC( J) )U1 ,U2,L1 
LI  PRINT  211»JfXC(J)  _.-  ,.» 

11  FCR.VAK  10X3HXG(  ,I2,LH)  =  ,  c  1  <; .  5  ) 

K  =  K  +  1 

i  o  r  omt  T  vi  I  C 


28 
29 

30 
31 
32 

33 
3L 
35 


36 

37 
38 


39 

LO 

21 C 


L2  CONTINUE      , 

IF(K)10U2,10L1,1CL2 

10L1  PRINT  2L0 
PRINT  2  02 
REAC  102, ITEST  _ 
ICHECK!  =  8HFCLC  TIM 
ICHECK2  -  8HPEAC  TIM 
IFCICHECKl  -  ITcST)U3,L5 ,L3 
IFUCHECK2  -  ITEST)LL,L5,LL 
PRINT  2CL 
STOP  8 
L5  READ  10L,TC,TF 


1CL2 


L3 
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FCRMAT(  10F8 

RLAD102, ITEST 

ICHECKl  =  8HCCMPUTE 

ICHECK2  =  8HFOLD  STE 
trucri/7  —  qudcap  ctc 


104  FCRMAK  10F8.U) 
46  RLA0102, ITEST 


51 
52 


INDIC1  =  0 

READ  104, (C( J) , J=1C3,1C9 ) 
IF(C(103)  -  TC  )52,53,52 
PRINT  204 


STOP  10 

53  IF(C( 104) )54,52,54 

54  NOT  =  1 

IF(C( 105)  -  TF)55,59,55 

55  IF(C( 106) )56,52,56 

56  NDT=  2 


56 

57 
58 


NDT=  2 

IF(C(107)  -TF)57,59,57 
IF(C(108) )58,52,58 
NDT  =  3 

IF(C(109)  -TF)52,59,52 
59  PRINT  212,TG,TF 
212  FORMAT131H  THE  TIMING  DATA  ARE  AS  FOLLOWS  ,// 
1         5X,  15HINITIAL  TIME  =  ,E11.5,/ 
5X,  15FFINAL  TINT 
IC1  -  1)61,60,61 
7  1  3 .  K  P 


2        5X,  15FFINAL  TIME    =  ,E1U5) 
IF( INDT 


It-  \   1  I\U1L  !   -   I  J  0  I  ,OU  ,  O  i 

60  PRINT  213, KP 
213  FORMAT(5X,35hTHE  STEP  SIZE  IS 

1         5X,31hSMALLEST  VARIABL 
AP  =  10.0**KP*1.0E-04/(TF  -  T 
DO  1060  J=1,N 
1060  A{J)  =  AP 

A( 100)  =  TF  -TC 

C(  104)  =  A{ lC0)»l.CE-05 

GO  TO  63 

61  NDTT  =  NDT«2 


NDTT  =  NDT*2 

DO  62  K=l , NDTT, 2 

PRINT  214,C{1C3+K),C(102+K),C( 10  4+K) 

FORMAT(5X,  15HSTEP  SIZE     =  ,E1 1 .5* 1  OH 

El  1.5) 
PRTNT  ?0? 


COMPUTED,  BASED  ON  ,/ 
E  OF  CRCER  l.OE  , 12) 
C) 


214  FORMATTsxIiSHSTEp'sIZE"   '^=' 7e1 lZsVlOH  FROM  T  =  ,E11.5,8H  TO  T  = 
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72 
73 


74 

75 
76 
77 


78 
79 
80 


81 
215 

82 
216 

83 
217 


1 


218 


1084 

1085 

1217 

84 

219 

85 
220 

86 
221 

241 

1087 
242 

1088 
243 

1089 
244 

1090 

245 

87 

2087 

2088 
2089 


223 


READ 
REAC 
ICHE 
I  F  (  I 
READ 
ICHE 
IF(  I 
ICHE 
I  F  (  I 
PR  IN 
STOP 
READ 
ICHE 
IF(  I 
IF(N 
IF(N 
NG2 
READ 
READ 
ICHE 
IF(  I 
PR  IK' 
FORM 
IF(N 
PRIN 
FORM 
GO  T 
PRIN 
FORM 

PRIN 
FORM 
DO  1 
IFt  I 
CONT 
GO  T 
PRIN 
FORM 
PRIN 
PRIN 
FORM 
IF(N 
PRIN 
FORM 
GO  T 
PRIN 
FORM 
PRIN 
FORM 
IF(N 
PRIN 
FORM 
IF(N 
PRIN 
FORM 
IF(N 
PRIN 
FORM 
IF(N 
PRIN 
FORM 
DO  2 
I  F  (  I 
CONT 
GO  T 
PRIN 
PRIN 
PRIN 
PRIN 
FORM 


102,  ( 
101,1 

CK  =  8 
CHECK 

102,  I 
CK  =  8 
CHECK 
CK  =  8 
CHECK 
T  204 
12 

1  0  1  ,  N 
CK  =  8 
CHECK 
G)79,8 
G  -  5) 
=  NG* 
1  104, 
101,  I 
CK  =  8 
CHECK 
T  215 
AT( 17H 
P)83,8 
T  216 
AT( 5X, 
0  84 
T  217, 
AT(5X, 
5X, 
T  218, 
AT{ 10X 
C84  J  = 
P(  J)  )  1 
INUE 
0  84 
T  1217 
AT  (/, 
T  202 
T  219 
AT(  17H 
G)86,8 
T  220 
AT(5X, 
0  2089 
T  221  , 
AT( 5X, 
T  241, 
AT  (10 
G  -  2) 
T  2  42, 
AT  (10 
G  -  3) 
T  243, 
AT  (10 
G  -  4) 
T  244, 
AT  (10 
G  -  5) 
T  245, 
AT  (10 
087  J  = 
G( J)  )2 
INUE 
0  2089 
1217 
200 
205, 
223, 


KT ITLE( J)„J=K,KP9) 
NCPR, ITEST 
H  INCREME 

-  ITEST)65,74, 65 
TEST 

HHCLD  GRA 

-  ITEST)75,81, 75 
HREAD  GRA 

-  ITEST)76,77, 76 


G,  ITEST 
H  GRAPHS 

-  ITEST)76,78, 76 
1,79 

80,80,76 

2 

(  IC(K  )  ,K=1 ,NG2 ) 

NCGR, I  TEST 

H  INCREME 

-  ITEST)76,81,76 

PRINT  SUMMARY ,/) 

2,63 

11FN0  PRINTOUT) 

INCPR 

I2,29H  INCREMENTS  BETWEEN  PRINTOUTS,/ 

25HTHE  VARIABLES  PRINTEC  ARE,/  ) 
(  IP( J  )  ,J  =  1  ,NP) 
,2FX( ,  13, 1H) ) 
1  ,NP 
084, 1C85, 1C84 


5X,22H*  X(0)  REPRESENTS  TIME  ) 


GRAPH  SUMMARY ,/) 

5,66 

9FN0  GRAPHS) 

INCGR 
I2,26H    INCREMENTS    BETWEEN    POINTS,/     ) 
I G  (  1  )  ,  I G  ( 2  ) 

X,13HGRAPH    A    IS    X(,I3,8H)     VS.    X(,I3,1H) 
87,  1087, 1C87 
I G  (  3  )  ,  I  G  (  4  ) 
X,13HGRAPH    B    IS    X(,I3,8H)     VS.    X(,I3,lH) 

87.1088. 1088 
I  G  (  5  )  ,  IG(6  ) 

X,13HGRAPH    C    IS    X(,I3,8H)     VS.    X(,I3,1H) 

87. 1089. 1089 
IG(7)  , IG ( 8  ) 

X,13HGRAPH    C    IS    X(,I3,8H)     VS.    X(,I3,lH) 

87.1090. 1090 
IG(9) , IG( 10) 

X,13HGRAPH  E  IS  X(,I3,8H)  VS.  X(,I3,1H) 

1  ,NG2 

087, 2C88, 2087 


(  ITITLEt J),J  =  1  ,5) 
NRC 
AT(5X,10FRUN  NUMBER,  I  2 , / // ) 
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IF(NP)88,90,88 

88  IF (NT) 89, 90, 89 

89  PRINT  224, (KTITlEf J) ,J=1 ,KP9) 

224  F0RMAT(2X,9( A8, ^X) ,A3) 
PRINT  201 

90  T  =  TO 

DT  =  C(  104) 
DO  9  1  J=l ,N 

91  X( J )  *  xO( j  ) 
LINES  =  0 
NCPTS  =  0 
NUMPTS  =  0 

301  IF(NP)302,312,302 

302  IF(N0PTS)303,308,303 

303  IF(XMCDF(NCPTS,bO*INCPR) )304,306,304 

304  IF( XMCDF( NOPTS,  10* INC PR)  )  30 5, 307,30  5 

305  IF(XMODF(NOPTS,    INCPR)  )312,308,312 

306  PRINT  200 

307  PRINT  2C1 

308  CALL  DERIV(T,X, XDOT,C) 
LINES  =  LINES  +  1 

DO  31 1  J=l ,NP 

IF(  IP( J) )310,309,310 

309  PR(J)  =  T 
GO  TO  311 

310  IPJ  =  IP( J) 
PR( J)  =  X{ IPJ) 

311  CONTINUE 

PRINT  225, {PR( J ) , J=l ,NP) 

225  FORMAT!  1 0 (  1  X  ,  E 1 1 . 5 ) ) 

312  IF(NG)313,318,313 

313  IF( XMCDF (NOPTS,  I  NCGR )  ) 31 8 , 3 1 4 , 3 1 8 

314  IF(XMODF(NOPTS,INCPR) ) 13 15, 1314, 1315 
1315  CALL  DERIV(T,X,XDOT,C) 

1314  DO  317  J=l ,NG2 

IF( IG( J) )316,315,316 

315  GRU)  =  T 
GO  TO  317 

316  IGJ  =  IG(J) 
GRU)  =  X(  IGJ) 

317  CONTINUE 

NUMPTS  =  NUMPTS  +  1 
Yl (NUMPTS)  =  GR( 1 ) 
XI  (NUMPTS)  =  GRU) 
Y2  (NUMPTS)  =  GRU) 
X2(NUMPTS)=  GR(4) 
Y3(NUMPTS)  =  GRU) 
X3(NUMPTS)  =  GR(6) 
Y4(NUMPTS)  =  GR(7) 
X4(NUMPTS)  =  GRU) 
Y5(NUMPTS)  =  GR(9) 
X5(NUMPTS)  =  GR( 10) 

318  NOPTS  =  NOPTS  +  1 

IFtLINES  -  250)1319,1318,1318 

1318  PRINT  1216 

1216  FORMAT  (//24F  STOP  AT  250  PRINT  LINES) 
GO  TO  341 

1319  IF(T  -  l.E  +  04)320,319,31  9 

319  PRINT  226 

226  FCRMAT(//19H  STCP  AT  T  =  10,000) 
GO  TO  341 

320  IFUQPTS  -  1C000  )322,321,321 

321  PRINT  227 

227  FCRMATI//26H    STCP    AT     10,000     INCREMENTS) 
GO    TO    341 

322  IFUUMPTS    -    900)324,323,323 

323  PRINT    228 

228  FCRMAT(//25H  STCP  AT  900  GRAPH  PCINTS) 
GO  TO  341 

324  IF(T  -  TF)326,325,325 
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229 
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337 
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340 
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341 

1341 
342 
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345 

346 

347 

348 

349 

350 
351 


PRINT  2 
FORMAT ( 
GO  TO  3 
DO  320 
IF( ABSF 
PRINT  2 
FORMAT* 
GO  TO  3 
CCNTINU 
IF( INOI 
STOP  13 
CALL  RK 
IF(NDT 
OT  =  C( 
GO  TO  3 
IFtNDT 
IF(T  - 
=  C( 
TO  3 


DT 

GO 


1 


352 


1 


353 


1 


354 


1 


355 


IF(T  - 
IF  (  T  - 
OT  =  C( 
T  =  T  +  D 
GO  TO  3 
CALL  RK 
GO  TO  3 
IF(NG) 1 
GO  TO  { 
ITITLE( 
GO  TO  3 
ITITLEt 
GO  TO  3 
ITITLEt 
GO  TO  3 
ITITLEt 
GO  TO  3 
ITITLEt 
GO  TO  3 
ITITLEt 
GO  TO  3 
ITITLEt 
GO  TO  3 
ITITLEt 
GO  TO  3 
ITITLEt 
LABEL  = 
MODCURV 
SFX  =  0 
SFY  =  0 
MINCFFX 
MINOFFY 
LABELNO 
MODE  = 
ITITLEt 
CALL  GR 

IFtNG  - 
ITITLEt 
CALL  GR 

IFtNG  - 
ITITLEt 
CALL  GR 

IFtNG  - 
ITITLEt 
CALL  GR 

IFtNG  - 
ITITLE! 


29 

//26H  NORMAL  STOP  AT  FINAL  TIME) 

4  1 

J=l ,150 

(X( J) )  -  l.E+04)328,327,327 

30,  J 

//11H  STCP  AT  X(,I3,1CH)  =  1C,000) 

4  1 

E 

CI  )331,332» 1331 

UTTA  (N,T,DT,X,C) 

-  1)334,333,334 
104) 

40 

-  2)337,335,337 

C(  105)  )333, 336,336 

106) 

40 

C( 105))333,333,338 

Ct 107) )336,339,339 

108) 

T 

01 

UTTA2  (N,T,CT,X,A  ,C) 

01 

341,356,1341 ■ 

342,314  3,344,345,3  46,347,348,349,35  0)  ,NRC 
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51 
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51 
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51 
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51 
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51 
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51 
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6) 
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8H 

8H 
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GH 
EH 
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RUN 
RUN 
RUN 
RUN 
RUN 
RUN 
RUN  7 
RUN  8 
RUN  9 
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=  1 
0 
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APH2 

2)3 
7)  = 
APH2 

3)3 

7)  = 
APH2 

4)3 
7)  = 

APH2 
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1 

SH  GRAPH  A 

(NU:VPTS,X1  ,Y  1,8,MCDCURV,LABEL,ITITLE,SFX,SFY, 
MINCFFX,  MINOFFY,  LABELNCMGDE) 
56,352,352 
EH  GRAPH  B 

(NUMPTS,X2,Y2,8,M0DCURV, LABEL, ITITLE, SFX, SFY, 
MINOFFX, MINOFFY, LA BELNC, MODE) 
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SH  GRAPH  C 
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EH  GRAPH  E 
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CALL  GRAPH2 


PRI 

IF( 

RET 

END 

SUB 

DIM 

CT( 

CT( 

CT( 

CT( 

DO 

TC 

DG 

XC( 

CAL 

DO 

AK( 

DO 

X(  J 

RET 

END 

SUB 


(NUMPTS,  X5,Y 
MINOFFX,MIN 

NT  200 

NRC  -  NR)300,357,357 

URN 
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R 
X( 
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0 

=  c. 
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=  1. 

1  =  1,  u 
T  +  CT{ 

J=l  ,N 

=  X(J) 
DERIV  ( 
J=l  ,N 
J)  =  DT 
J=1,N 
=  X(  J  ) 


KUTTA  (N,T,D 
150)  ,AKU,  10 


I)*DT 

+  CT(  I) *AK( 
TC,XC,XDOT,C 


5,3,M0DCURV,LABEL, ITITLE,SFX,SFY, 
CFFY,LAEELNC,MODE) 
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I-l.J) 
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DIME 
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DO  1 
X2(  I 
X22( 
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XS(  I 
TS  = 
CALL 
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DO  6 
E21 
E21R 
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IF 

IF 
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GO 
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URN 
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RKUTTA 
DT 

RKUTTA 
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RKUTTA 
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THE 
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OF 
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METHOD 
-360).  I 
F  THE  IN 
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TA,  WHIC 
GUEMENTS 

EGUATIO 
ATICN  ST 
^PLETION 
ACH  STEP 
ES  (ALSO 

ERROR  P 
LES. 

AL  TIME. 
DUCED  FR 
IFIED  AC 


CESCRIBEC  IN  COMM.  OF 
T  CHECKS  THE 
TEGRATICN  AND  ADJUSTS 
UNGE-KUTTA  INTEGRATION 
H  IS  AVAILABLE  IN  THE 

ARE, 
NS.  MAXIMUM  N  =99. 
EP  (UPDATED  BY  THIS 

CF  EACH  STEP  )  . 

(ALSC  UPDATED). 

UPDATED)  . 
ER  UNIT  TIME  FOR  EACH 


CM  TFE  VALUE  STATED 
CURACY  HAS  BEEN 


, XS(  100)  ,X22(  1C0),X2( 100),A(  ICC)  ,C( 150) 


(N,TS,H,X22, C) 
(N,TS,H,X2,C ) 

) 


(N,TS,H,XS,C ) 


CONDITI 
E21R  /  ( 
U  -  U2)6 

U 
)  =  XS( I 
U2  -  1.0 
H  -  A( 10 

A (100)* 

SETS  TH 
C  16 
0  1=  1  ,N 


)  -XS( I ) )*0. 
ABSF(E21 ),AB 
-1 
CN  PREVENTS 
A(  I)*6.0*H) 
,6,5 

)  -  E21 

)11 ,7,7 

C)*l  .0E-C9)8 

l.CE-09 

E  MINIMUM  ST 


C6666666667 

SF( ( ABSF(XS( I ) )+ABSF(X22( I ) ) 

.99*ABSF(X( I  )  )  )*1  .0E-C9)  ) 

RCUNC-OFF  ERROR  FROM  TAKING  CONTROL 


,8,9 

EP  SIZE  TO  1.0E-09  TIMES  THE  TCTAL  TIME. 
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X22(  I)  =  X2(  I) 

10  X2( I )  =  X( I ) 
TS  =  T 

H  =  C.5*H 

GO  TO  2 

THIS  RECYCLES  THE  INTEGRATION  IF  THE  TRUNCATION  ERROR  IS  EXCESSIVE 

11  IF  (U2  -  0.031  ) 12,  13,13 

12  DT  =  2.0*H 
GO  TO  14 

13  DT  =  SQRTF( SCRTF ( 0 . 5/U2)  )*H 

14  IF  (DT  -  A(1C0)*0. 001 ) 16  ,16,15 

15  DT  =  A( 100)*C.CC1 

THIS  SETS  THE  MAXIMUM  STEP  SIZE  TO  1.0E-03  TIMES  THE  TOTAL  TIME. 

16  T  =  T  +  2.0*1- 
C(  104)  =  DT 
DO  17  I  =1,N 

17  X( I )  =  XS{ I ) 

THIS  UPDATES  T  AND  X(I).  DT  IS  UPDATED  BY  STATEMENT  8,12,13  OR  15. 

RETURN 

END 

FUNCTION  RELAY  (R,DZONE,V) 

IF( ABSF (R)  -  DZCNE)3,3,4 

3  RELAY  =  0.0 
RETURN 

4  RELAY  =  SIGNF(V,R) 
RETURN 

END 

SUBROUTINE  OERIV  (T,X,XDCT,C) 

DIMENSION  Xt  150)  ,XDOT( 10  0),C( 150) 
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COMMENTS 

C      SUBROUTINE  FCR  OPEN  -  CLCSED  CYCLE  SYSTEM  USING  A  TIMING  CKT  WHERE 
C      TIME  TO  SWITCH  AND  TIME  TO  TURN  OFF  IS  PROP  TC  THE  COMMAND  SIGNAL 
C      AND  THE  RELAY  OUTPUT  IS  A  CONSTANT  IN  THE  NONLINEAR  REGION 

R  =  C  (  1  ) 

ERROR  =  R  -  X( 1 ) 

IF(SWl)  599,  599,  6CC 

599  IF  (ABSF(ERRCR  -  C(2)))  603,  60C,  6C0 

600  SIGN  =  ABSF( ERRCR)/ERROR 
SWl  =  1.0 
SW2  =  O.C 
IF  (T  -  TS  -  TIME)  601,  602,  602 

601  DRIVE  =  C(6)  *  SW1 
GO  TO  604 

602  DRIVE  =  -C(6)*  SWl 
IF(ABSF(XDCT( 1 ) )  -  .05)'  603,  603,  604 

603  SWl  =  0.0 
SW2  =  1.0 
TIME  =  X(0) 

604  POWER  =  C(3)  *  ERROR  *  SW2  +  DRIVE  *  SWl 
XDOTM  )  =  X(2) 

XD0T(2)  =  POkER  *  C(4)  -  X(2)  *  C(5) 
TS  =  0.075  *  ABSF(R)  +  0 .05 
TO  =  0.092  *  ABSF(R)  +  0.105 
X(3)  =  SWl 
X(4)  =  SW2 
X(5)  =  DRIVE 
X(6)  =  POWER 
X<  101  )  =  ERRCR 
X(  102)  =  -XDCT( 1 ) 
RETURN 

COMMENTS  THIS  IS  BLUE  DECK  ONE 

C      X{1)=SYSTEM  CUTPUT,  X(2)=SYSTEM  RATE,  C(l)=  INPUT,  C(2)  =1/2WIDTH 
C      OF  LINEAR  ZONE,  C(3)=SYSTEM  GAIN,  C(4)=PLANT  GAIN,  C(5)=1/TIME  C 
C      C(6)=MAGNITUCE  CF  DRIVING  VOLTAGE  IN  NONLIN  REGION 

END 

END 


Fig*    4l«      Data   input  for   Program  Analog. 


KIlltKMAN     f^ALLtY      IHtbIS 
02    EQUATIONS 
04    RUNS 

01     TITLE    CARDS 

OPEN-CLOSED    LOOP    SERVC    SYSTEM    WITH    SWITCHING    =    F(TIME,     INPUT) 

ZERO    COEFFS 

06  COEFF.  CARCS 


1 

0.5 

2 

0.25 

3 

40.00 

,  4 

1  .0 

5 

5.0 

6 

100.0 

ZERO 

ICS 

00  IC 

CARDS 

READ 

TIME  DATA 

0.0 

1.5 

COMPUTE  STEP  SIZE 

-1  IS  ORDER  OF  SMALLEST  VARIABLE 

READ  PRINT  DATA 

07  VARIABLES  PRINTED 

0  12    3    4    5    6    7 
02  TITLE  CARDS 

TIME   SYSTEM   SYSTEM   SWITCH   SWITCH    CRIVE    PCWER 
OUTPUT    RATE      ONE      TWO 
20  INCREMENTS  BETWEEN  PRINTS 
READ  GRAPH  DATA 
02  GRAPHS 

1  0  102  101 

01  INCREMENTS  BETWEEN  POINTS 

HOLD  COEFFS 

01  COEFF.  CARDS 

1   1.0 
ZERO  ICS 

00  IC  CARDS 
HOLD  TIME  DATA 
HOLD  STEP  SIZE 
HOLD  PRINT  DATA 
HOLD  GRAPH  DATA 
HOLD  COEFFS 

01  COEFF.  CARDS 
1   1.5 

ZERO  ICS 

00  IC  CARDS 
HOLD  TIME  DATA 
HOLD  STEP  SIZE 
HOLD  PRINT  DATA 
HOLD  GRAPH  DATA 
HOLD  COEFFS 

01  COEFF.  CARCS 
1   2.0 

ZERO  ICS 
00  IC  CARDS 
HOLD  TIME  DATA 
HOLD  STEP  SIZE 
HOLD  PRINT  DATA 
HOLD  GRAPH  DATA 


=  72 


Fig.  42.   Typical  outnut  for  Program  Analog. 


OPEN-CLOSED  LOOP  SERVO  SYSTEM  WITH  SWITCHING  =  FtTIME,  INPUT) 
k    RUNS  ARE  CALLED  FOR. 


RUN  NUMBER  1 
THE 


NON-ZERO  DATA 

COEFFICIENTS  ARE  AS  FOLLOWS 

C(  1)  = 

.50000E+03 

C(  2)  = 

.25000E+00 

C(  3)  = 

.UOOOOE+02 

C(  4)  = 

. 10000E+01 

C(  5)  = 

.50000E+01 

C(  6)  = 

.I0000E+03 

THE  NON-ZERO  INITITAL  CONDITIONS  ARE  AS  FOLLOWS 
NONE 

THE  TIMING  DATA  ARE  AS  FOLLOWS 

INITIAL  TIME  =   .OOOOOE+OO 

FINAL  TIME    =   .15000E+D1 

THE  STEP  SIZE  IS  COMPUTED,  BASED  ON 

SMALLEST  VARIABLE  OF  ORDER  l.OE-1 

PRINT  SUMMARY 

5  INCREMENTS  BETWEEN  PRINTOUTS 
THE  VARIABLES  PRINTED  ARE 

X(  0) 

X(  1  ) 

X<  2) 

X(  3) 

X(  4) 

X(  5) 

X(  6) 

*  X(0)  REPRESENTS  TIME 

GRAPH  SUMMARY 

NO  GRAPHS 
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Fig.  42.  Oont. 


KITTERMAM 
RUN  NUMBER 


MALLEY 
1 


THESIS 


TIME 


SYSTEM 

OUTPUT 


SYSTEM 

RATE 


SWITCH 

ONE 


SWITCH 
TWO 


DRIVE 


.0000 
.9300 
.1281 
.2781 
.U281 
.5781 
.7281 
.8781 
.1028 
.1178 


OE  +  00 
OE-03 
0E-01 
OE-01 
OE-01 
OE-01 
OE-01 
OE-01 
1E  +  00 
1E  +  00 


OOOOOE+OO 
43173E-04 
80324E-02 
36938E-01 
85432E-01 
15210E+00 
23563E+00 
33479E+00 
44847E+00 
57562E+00 


.OOOOOE+OO  .OOOOOE+OO  .10000E+01  -. 10000E+03 

.92784E-01  .10000E+01  .OOOOOE  +  OO  -.  1  0000E+03 

.12408E+01  .10000E+01  .OOOOOE+OO  -.  10000E  +  03 

.25963E+01  .10000E+01  .OOOOOE+OO  -. 10000E+03 

.38538E+01  .10000E+01  . OOOOOE+CO  -. 1 0000E+03 

.50205E+01  .10000E  +  01  .OOOOOE -.  30  -.  1  0000E  +  03 

.61029E+01  .10000E+01  .OOOOOEvOO  -.10000E+03 

.71370E+01  .10000E+01  .OOOOOE+OO  -.10000E+03 

.80386E+01  .10000E+01  . OOOOOE +00  -, 1 0000E+03 

.89329E+01  .10000E+01  .OOOOOE+OO  -. 1 0000E+03 


13281E+00 
14500E+00 
14500E+00 
14500E+00 
14500E+00 
14502E+00 
14579E+00 
15715E+00 
17215E+00 
•18715E  +  00 


.2021 
.2171 
.2321 
.2U71 
.2621 
.2771 
.2921 
.3071 
.3221 
.3371 


.3521 
.3671 
.3821 
.3971 
.U121 
.4271 
.4421 
.4571 
.4721 
.4871 


5E  +  00 
5E  +  00 
5E  +  00 
5E  +  00 
5E  +  0C 
5E  +  00 
5E  +  00 
5E  +  00 
5E  +  00 
5E  +  0C 


5E  +  00 
5E  +  00 
5E  +  00 
5E+00 
5E  +  00 
5E  +  00 
5E  +  00 
5E  +  00 
5E  +  00 
5E  +  00 


.5021 
.5171 
.5321 
.5471 
.5621 
.5771 
.5921 
.6071 
.6221 
.6371 


5E  +  00 
5E  +  00 
5E  +  00 
5E  +  00 
5E  +  00 
5E  +  00 
5E  +  00 
5E  +  00 
5E  +  00 
5E  +  00 


.71525E+00 
.83727E+00 
.83730E+00 
•.83730E  +  00 
•.83731E  +  00 
■.83755E  +  00 
•.84538E  +  00 
•.95485E  +  00 
•.10790E  +  01 
•.1  1805E  +  01 


.12601E+01 

.  13187E+01 

.13574E+01 

.13771E+01 

. 1 3791 E+01 

•.13648E+01 

.13353E+01 

-.12920E+01 

•.12363E+01 

-.11696E+01 


.10933E+01 

.10087E+01 
.91707E+00 
■.81983E  +  00 
.71817E+00 
.61331E+00 
■.50638E  +  00 
■.39845E  +  00 
■.29054E+00 
■.18358E  +  00 


.78422E-01 

.24149E-01 
. 12343E+00 
.21880E+00 
.30972E+00 
.39571E+00 
.47640E+00 
.55146E+00 
.62065E+00 
.68380E+00 


.97D 
.  103 
.103 
.  103 
.103 
.103 
.  102 
.904 
.751 
.602 


48E+01 
13E+02 
14E+02 
14E+02 
13E+02 
11E+02 
31E+C2 
69E+01 
48E+01 
82E+01 


.459 
.322 
.193 
.712 
.425 
.147 
.244 
.331 
.409 
.473 


69E+01 
91E+01 
22 E+01 
32E+00 
63E+00 
77E+01 
09E+01 
34E+01 
41E+01 
29E+01 


.53305E+01 
.58334E+01 
.63385E+01 
.66436E+01 
.68969E+01 
•70722E+01 
.71735E+01 
.72053E+01 
.71725E+01 
.70798E+01 


.69325E+01 
.67358E+01 
.64949E+01 
.621 52E+01 
.593 18E+01 
.55599E+01 
.51946E+01 
.48107E+01 
.441 28E+01 
.40354E+01 


.10000E+01 
.10000E+01 
.100C0E+01 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 


.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 


E  +  00 
E  +  00 
E+00 
E  +  00 
E  +  00 
E  +  00 
E  +  00 
E  +  00 
E  +  00 
E  +  00 


.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 


.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
•.OOOOOE  +  OO 
.OOOOOE+OO 
•OOOOOE+OO 
.OOOOOE+OO 


.0000 
.0000 
.0000 
.1000 
.  1000 
.  1000 
.1000 
.1000 
.1000 
.  1000 


OE  +  OO 
OE  +  00 
OE  +  OO 
OE  +  0" 
OE+0 
OE  +  O 
OE  +  0 
OE+O 
Ot  +  0 
OE  +  O 


10000E+0 
10000E+0 
10000E+0 
10000E+0 
10000E+0 
10000E+0 
10000E+0 
10000E+0 
10000E+0 
10000E+0 


10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 


E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E+0 
E  +  0 
E  +  0 
E  +  0 


.10000E+0 
.10000E+0 
.10000E+0 
.10000E+0 
. 10000E+0 
.10000E+0 
.10000E+0 
. 10000E+0 
.10000E+0 
.10000E+0 


10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 


E  +  03 
E  +  03 
E  +  03 
E  +  03 
E  +  03 
E  +  03 
E  +  03 
E+03 
E+03 
E  +  03 


10000E+03 
10000E+03 
10000E+03 
10000E+03 
10000E+03 
10000E+03 
10000E+03 
10000E+03 
10000E+03 
10000E+03 


10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 


E  +  03 
E  +  03 
E+Cr 
E+Oo 
E  +  03 
E  +  03 
E  +  03 
E+0  3 
E  +  03 
E+03 


10000E+03 
10000E+03 
10000E+03 
10000E+03 
10000E+03 
10000E+03 
10000E+03 
10000E+03 
10000E+03 
10000E+03 


I    EX  response,  second  «***- 
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